Proteomic insight into human directed selection of the domesticated chicken Gallus gallus.

Chicken domestication began at least 3,500 years ago for purposes of divination, cockfighting, and food. Prior to industrial scale chicken production, domestication selected larger birds with increased egg production. In the mid-20th century companies began intensive selection with the broiler (meat) industry focusing on improved feed conversion, rapid growth, and breast muscle yield. Here we present proteomic analysis comparing the modern broiler line, Ross 708, with the UIUC legacy line which is not selected for growth traits. Breast muscle proteome analysis identifies cellular processes that have responded to human directed artificial selection. Mass spectrometry was used to identify protein level differences in the breast muscle of 6-day old chicks from Modern and Legacy lines. Our results indicate elevated levels of stress proteins, ribosomal proteins and proteins that participate in the innate immune pathway in the Modern chickens. Furthermore, the comparative analyses indicated expression differences for proteins involved in multiple biochemical pathways. In particular, the Modern line had elevated levels of proteins affecting the pentose phosphate pathway, TCA cycle and fatty acid oxidation while proteins involved in the first phase of glycolysis were reduced compared to the Legacy line. These analyses provide hypotheses linking the morphometric changes driven by human directed selection to biochemical pathways. These results also have implications for the poultry industry, specifically Wooden Breast disease which is linked to rapid breast muscle growth.


Introduction
The chicken has played an important role in human culture and nutrition since its domestication [1,2]. The majority of the modern chicken's genome was derived from the red junglefowl with documented introgression of alleles from the grey, Ceylon, and green junglefowl [3]. Early domestication likely led to larger birds that matured quicker and increased egg production in comparison with their wild progenitors [4,5]. The advent of modern agriculture in the early 20 th century led to intensive genetic selection for either meat (broiler) or egg (layer) production traits in the chicken. Initial efforts selecting for increased broiler mass led to adipose tissue accumulation, probably because the selection did not channel the increased level of nutrients to a particular tissue [6]. Ultimately, improving broiler production traits led to selection for a combination of larger breast muscle mass, improved feed efficiency and rapid growth.
Comparative studies describing differences between selected and unselected (legacy) lines provide insight into the impact of artificial selection of species. Several studies have characterized the differences between modern broilers and legacy lines that have not been subjected to production level human directed selection [7][8][9][10][11][12][13]. In our work [14][15][16] we compared the modern Ross 708 broiler line and the legacy University of Illinois, Urbana Campus (UIUC) lines (Fig 1) [17,18]. In the legacy line the breast muscle comprises approximately 9% of the body mass, while in the modern broiler this tissue constitutes up to 22% of the body mass [14]. Evolving in the tropics, the red junglefowl had little need for long distance flying. The wild chicken is an episodic flier, only needing the ability to fly up into a tree to escape predators and to roost. Like the legacy line, the breast muscle of the red junglefowl constitutes approximately 9% of its body mass [19]. For comparison, the breast muscle of birds capable of more sustained flight averages 17% while that of hummingbirds varies between 25-30% [20]. The increase in the modern broiler's breast muscle mass, which almost reaches the hummingbird level, came at the expense of other tissues. For example, the normalized masses of the heart, spleen and brain are larger in legacy lines compared with modern broilers [13,14,21]. In modern broilers, this likely is responsible for increased incidence of cardiomyopathy [22,23] and immune deficiencies [10]. The reduced brain mass in broilers could contribute to the behavioral differences seen in modern broilers compared with other lines. These include a reduction in movement, along with reduced fear and risk aversion in comparison with legacy or layer lines [24]. In addition, a myopathy called Wooden Breast disease has arisen during the selection process [25][26][27]. This myopathy is apparently due to selection for rapid growth and is characterized by necrosis, fibrosis, and immune cells infiltration that results in a product that is unappealing to consumers [28][29][30][31]. This disease results in significant losses to the poultry industry as the meat is condemned at processing [32].
Our previous research involved comparing the breast muscle transcriptome of 6-day old (D6) modern and legacy lines [16]. That comparison identified differently expressed genes affecting growth factors, lipid metabolism and the pentose phosphate pathway. In particular, the modern transcriptome indicated elevated expression of insulin like growth factor 1 (IGF1) which stimulates muscle hypertrophy combined with reduced expression of myostatin (MSTN), an inhibitor of skeletal muscle hyperplasia [33]. These studies suggest that human directed selection increased growth factor stimulation in the modern breast muscle. Furthermore, elevated expression of genes involved in the pentose phosphate pathway and lipid metabolism would provide necessary intermediates and energy required for rapid growth.
Here we report a comparative study between modern Ross 708 and legacy UIUC breast muscle proteomes. The Ross 708 modern line is a commercial breed that is produced for human consumption. The breeding process for this line is proprietary information owned by Aviagen. As a result, a detailed description of the breeding process and the specific genetic background of Ross 708 is not publicly available. The UIUC legacy line birds were obtained from the University of Illinois, Urbana Campus. They are the progeny of a cross between males from the inbred New Hampshire and females from the inbred Columbian lines [14]. The results support and extend the conclusions of our prior studies, providing further insight into the impact of selection that yielded modern broilers [14,16]. The data also have implications for the negative effects of intense selection on the broiler chicken along with the emergence of Wooden Breast disease [25,28,29,34,35].

Animal care and sample collection
Animal raising, handling and sample collection methods were approved by the University of Delaware Institutional Animal Care and Use Committee (Permit Number: 2703-12-10). Six male UIUC and six male Ross 708 breast muscle samples were collected at day 6 (D6) posthatch, frozen in liquid nitrogen and stored at -80˚C until processed for proteome analysis. All samples were taken from the posterior region of the left pectoralis major muscle [16]. We chose D6 because this is the earliest day we detect difference in the normalized breast muscle mass between the two lines [14].

Proteomic analysis
For each of the muscle samples, six technical replicates were analyzed by mass spectrometry. In each case 50 mg of each muscle sample was subjected to differential detergent fractionation and 20 μg of each fraction was trypsin digested as previously described [36,37]. Following digestion, each fraction was desalted using a peptide macrotrap (Michrom BioResources) according to the manufacturer's instructions. After desalting, each fraction was further purified using a strong cation exchange macrotrap to remove any residual detergent, which could interfere with the mass spectrometry. Fractions were dried and resuspended in 10 μl of 2% acetonitrile, 0.1% formic acid and transferred to low retention vials in preparation for separation using reverse phase liquid chromatography.
An Ultimate 3000 (Dionex) high performance liquid chromatography system coupled with an LTQ Velos Pro (Thermo) mass spectrometer were used for peptide separation and mass spectrum acquisition. The U3000 was operated at a flow rate of 333 nl per minute and equipped with a 75 μm x 10 cm fused silica column packed with Halo C18 reverse phase material (Mac-Mod Analytical). Each peptide sample was separated using a 4 h gradient from 2% to 50% acetonitrile with 0.1% formic acid as a proton source. The column was located on a Nanospray Flex Ion Source (Thermo) and connected directly to a silica Nanospray emitter to minimize peak broadening. High voltage was applied using a stainless-steel junction between the column and the emitter. Scan parameters for the LTQ Velos Pro were one MS scan followed by 10 MS/MS scans of the 5 most intense peaks. MS/MS scans were performed in pairs, one using collision induced dissociation (CID) and the other using higher-energy collisional dissociation (HCD). Dynamic exclusion was enabled with a mass exclusion time of 3 min and a repeat count of 1 within 30 sec of initial m/z measurement.

Protein identification
Spectrum matching programs X!tandem [38] and OMSSA [39] were used via the University of Arizona High Throughput Computing Center. Raw spectra were converted to MGF format for analysis using the MSConvert tool from the ProteoWizard software suite [40]. X!tandem was run with 12 threads, precursor and fragment tolerance of 0.2 Da, and up to two missed tryptic cleavages. Variable modifications used in the searches were: single and double oxidation of Methionine, carbamidomethylation of Cysteine, and phosphorylation of Tyrosine, Threonine, and Serine. OMSSA was run with 12 threads, precursor and fragment tolerance of 0.2 Da, up to two missed tryptic cleavages, and set to XML output format. A custom Perl script was used to parse XML search results from both X!tandem and OMSSA. Peptides with evalues � 0.05 were accepted and single spectrum identifications were rejected unless they were identified by both search engines. To verify data set quality, decoy searches were performed in the exact manner as before, but with a randomized version of the protein databases. False discovery rates ranged from 0.9-2.3% with an average of 1.7%.

Identifying differentially expressed proteins
Differential expression of proteins between lines was performed pairwise using peptide elution profiles as described previously [41]. Precursor mass spectra were extracted from the raw data in MS1 format using the MSConvert software from the ProteoWizard toolset. Peptide precursor m/z values were extracted from the previously compiled protein identifications using Perl. Peptide intensities were summed for each protein on a per-replicate basis. Data were normalized based on the mode of each replicate rather than the mean to minimize the effect of extreme values. A resampling analysis was performed for each pairwise comparison. Proteins were considered to be differentially expressed if the difference in means between conditions resulted in a P-value � 0.05.

Functional analysis
Gene Ontology (GO) [

Data deposition
Transcript data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus (GEO) and are accessible through GEO Series accession number GSE65217. Proteomics data is available from ProteomeXchange (PXD005288).

Proteomic results
In the comparison between D6 modern Ross 708 and legacy UIUC breast muscle, a total of 222 differentially expressed proteins were detected (P-value <0.05), with 173 enriched in modern samples and 49 enriched in legacy samples. Among the 222 differentially expressed proteins, 130 (57%) exhibited the same direction of enrichment as shown in an earlier transcriptome study [16] (S1 Table).

Gene Ontology [GO) [42].
Muscle from legacy birds were enriched for GO terms relevant to myofibers and energy production in striated muscle, along with vesicle transport (Fig 2). Muscle from modern birds were enriched for terms including stress response, myofibers, translation, energy production, metabolism, vesicle transport and innate immunity.

Stress response
In the comparison between modern and legacy muscle, 49 differentially regulated proteins were part of the stress response. Five stress response proteins were elevated in the legacy breast muscle including cytochrome c oxidase subunit 6A1 (COX6A1), capping actin protein of muscle Z-line subunit beta (CAPZB), ribosomal protein S27 (RPS27), G protein subunit beta 4 (GNB4) and DnaJ heat shock protein family (Hsp40) member B11 (DNAJB11). 44 stress-associated proteins were elevated in the modern chicken muscle compared to the legacy samples. These function in a variety of stress related processes including acting as chaperones and cochaperones, responding to genotoxic stress, redox regulation, and protein degradation. Among these were four proteasome subunits that likely offset the increased burden of misfolded proteins due to the large amount of protein synthesis required for muscle hypertrophy. Nine heat shock proteins (HSPs) were also elevated in muscle from the modern line, all of which function as chaperones controlling the proper folding of client proteins. Five members of the TCP complex that acts as an ATP dependent chaperone controlling actin folding were enriched in the modern samples [5,6]. Four isomerases responsible for controlling disulfide or proline conformation and three proteins that control redox stress were also elevated in the modern line.

Translation
The most striking contrast between modern broilers such as Ross 708 chickens and earlier breeds is the difference in both normalized and total breast muscle yield. Some of the increase in normalized breast muscle of modern broilers can be attributed to hypertrophy due to increased protein synthesis [54,55]. This is consistent with the enrichment in the Ross 708 birds for protein translation initiation factors and ribosomal structural proteins (S1 Table and Fig 2). In addition to ribosomal proteins, five enzymes encoding tRNA ligases were enriched in Ross 708 breast muscle including: lysyl-tRNA synthetase 1 (KARS1), seryl-tRNA synthetase (SARS), glycyl-tRNA synthetase 1 (GARS), alanyl-tRNA synthetase (AARS) and glutamyl-prolyl-tRNA synthetase 1 (EPRS). These would support the elevated translation required for the muscle hypertrophy seen in the modern line.

Glycolysis
Muscle from the UIUC legacy line is enriched for two proteins involved in glycolysis: hexose kinase 1 (HK) and muscle phosphofructokinase (PFKM) (Fig 3). HK1 catalyzes the first reaction of glycolysis by phosphorylating glucose while PFKM drives the commitment step to glycolysis. Elevated levels of these glycolytic enzymes are consistent with the fast twitch nature of breast muscle from wild chickens, which is driven by glycolysis.

Pentose Phosphate Pathway (PPP)
The modern line breast muscle contained higher levels of the PPP enzymes transketolase (TKT) and phosphoribosyl pyrophosphate synthetase (PRPS1 and PRPS2). These enzymes function in the nonoxidative phase of PPP returning fructose 6 phosphate and glyceraldehyde 3 phosphate to glycolysis and providing ribose 5 phosphate as a precursor for purine and pyrimidine synthesis (Fig 3).

TCA cycle
The main pathway affecting energy production enriched in the modern line breast muscle was the TCA cycle (Fig 3). Most enzymes directly involved in the TCA cycle were elevated in the modern line muscle, including citrate synthase (CS), aconitase 2 (ACO2), isocitrate dehydrogenases (IDH, IDH3A, IDH3B), dihydrolipoamide S-succinyltransferase (DLST, a component of α-Ketoglutarate dehydrogenase), succinyl-CoA ligases (SUCLA2 and SUCLG1) and fumarate hydratase (FH). The components not elevated were succinate dehydrogenase complex flavoprotein subunit A (SDHA) and malate dehydrogenase (MDH). The enrichment for components of the TCA cycle indicates the central role of this pathway in the energy metabolism of the modern Ross 708 chickens compared with the legacy UIUC birds. In addition to providing coenzymes for oxidative phosphorylation, the TCA cycle provides intermediates for cataplerosis to be used as precursors for nucleotides, amino acids, and lipids. Anaplerosis could be supported by elevated levels of glutamic-oxaloacetic transaminase 2 (GOT2) that converts glutamate to α-ketoglutarate thus allowing for replenishing TCA cycle metabolites.

Beta-oxidation
Enzymes involved in fatty acid beta oxidation are enriched in the modern birds compared with the legacy line (Fig 4). One enriched enzyme is acyl-CoA dehydrogenase family member 9 (ACAD9), the rate limiting enzyme in the oxidation of fatty acyl CoA. ACAD9 is responsible for introducing trans double bonds into palmitoyl-CoA and initiating the beta-oxidation of this common lipid. Also enriched is enoyl-CoA delta isomerase 1 (ECI1) which is necessary for beta oxidation of unsaturated fatty acids. The transcript encoding hydroxyacyl-CoA dehydrogenase (HADH) is also elevated in the modern Ross 708 muscle. This enzyme acts repeatedly on lipids, sequentially removing two carbon units by oxidizing a 12-carbon fatty acid to acetoacetyl-CoA. Acetyl-CoA acyltransferase 2 (ACAT2), also enriched in Ross 708 breast muscle, oxidizes acetoacetyl-CoA to acetyl-CoA. The product of beta oxidation, acetyl-CoA, can then be metabolized via the TCA cycle to generate energy, or used in anabolic reactions to support rapid breast muscle growth.

Amino acid metabolism
In addition to TCA enzymes affecting amino acid anaplerosis, modern skeletal muscle was enriched for enzymes affecting glycine, lysine, and methionine metabolism. Enzymes affecting glycine impact choline and creatine levels include aldehyde dehydrogenase 7 family member A1 (ALDH7A1), betaine-homocysteine S-methyltransferase (BHMT), dimethylglycine dehydrogenase (DMGDH) and glycine amidinotransferase (GATM). ALDH7A1, BHMT and DMGDH form a pathway in the conversion of choline to sarcosine that is found in high concentrations in skeletal muscle. Choline is a precursor to phosphatidylcholine, an important component of cellular membranes. GATM is part of the pathway that mediates the interconversion between creatine and glycine. Creatine is also found at high concentration in muscle as it is important for energy storage as creatine phosphate. Further supporting creatine phosphorylation is creatine kinase B (CKB), which is elevated in modern breast muscle compared with the legacy line tissue. Three enzymes enriched in the modern line are involved in S-adenosylmethionine (SAM) production from methionine. These include methionine adenosyltransferase (MAT1A), aldehyde dehydrogenase 7 family member A1 (ALDH7A1), betainehomocysteine S-methyltransferase (BHMT), and adenosylhomocysteinase (AHCY). SAM functions as a universal methyl donor in biological systems.

Electron transport chain
Enzymes associated with the mitochondrial respirasome showed enrichment in the legacy line breast muscle. These enzymes include NADH:ubiquinone oxidoreductase core subunit S7 (NDUFS7), cytochrome c reductases (UQCRFS1, UQCRC1, UQCRC2), cytochrome c oxidase subunit 6A1 (COX6A1) and ATP synthase F1 subunit gamma (ATP5C1), a component of the ATP synthase central stalk. Four mitochondrial components were found enriched in modern Ross 708 muscle: one subunit of the electron-transfer-flavoprotein (ETFA) along with inorganic pyrophosphatase 2 (PPA2). One NADH dehydrogenase (NDUFS4) and one component of the ATPase stalk, ATP synthase peripheral stalk subunit d (ATP5H) were also elevated in the modern muscle.

Vesicle and protein transport
The legacy line breast muscle was enriched for proteins affecting vesicle transport including: COPI coat complex subunit beta 2 (COPB2), transmembrane p24 trafficking protein 2 (TMED2), and amphiphysin (AMPH). COPB2 functions in retrograde transport from the Golgi to the ER and is involved in recognition of specific vesicle cargo proteins [56]. The TMED proteins function in anterograde transport of vesicles from the ER to the Golgi apparatus and play roles in cargo selection [9]. AMPH is enriched in neural tissue where it is involved in endocytosis through interaction with clathrin [57,58]. Some instances of the paraneoplastic Stiff-person Syndrome (SPS) [59,60] are due to AMPH autoantibodies produced in breast and lung cancer patients. SPS is characterized by muscle spasms, rigidity and hypertrophy that arise from the effect of the autoantibody on nerve cells that control muscles [61]. Modern Ross 708 breast muscle was enriched for vesicle regulatory components including: cell division cycle 42 (CDC42), GDP dissociation inhibitor 2 (GDI2), Rac family small GTPase 1 (RAC1) and ras homolog family member A (RHOA).

Innate immunity
Gene Ontology analysis identified 39 proteins that function in innate immunity. Two, gelsolin (GSN) and bridging integrator 2 (BIN2) were enriched in the legacy line breast muscle. Elevated expression of GSN is associated with a decreased response to a variety of inflammatory stimuli [62][63][64][65][66][67] and elevated BIN2 levels are associated with decreased phagocytic activity [68]. The 26 proteins enriched in the modern Ross 708 breast muscle are associated with neutrophil degranulation. For example, cathepsin B CTSB is secreted during neutrophil degranulation and degrades collagen upon release into the extracellular space [69,70]. Peroxiredoxin 6 PRDX6 and RAC1 associate with and activate NADPH oxidase, a major source of reactive oxygen species [71][72][73]. RAC1 [74], RHOA [75] and WD repeat domain 1 (WDR1) [76] participate in the polarization that drives neutrophil migration. Birds do not have cells named neutrophils, but neutrophil function is carried out by heterophils in avian species [77]. It is reasonable to hypothesize that the innate immunity proteins are present in heterophils within the breast muscle.

Discussion
The differences between modern broilers (Ross 708) and legacy lines (UIUC), arose from human directed selection leading to broiler's rapid growth, improved feed efficiency and increased breast muscle mass. Rapid growth is evidenced in the time from hatch to market. Typically, for the legacy line it takes 16 weeks for birds to reach market weight, while birds from the modern line take seven weeks (see Fig 1 for growth comparison). The improved feed efficiency arose in part from selection lengthening the absorptive segments of the small intestine combined with earlier maturation of the liver [14]. Selection for increased breast muscle has generated modern birds with more than twice the breast muscle mass of legacy lines [14]. The proteomic results presented here indicate that selection caused metabolic reprogramming that supports the excessive breast muscle growth seen in the modern lines.
In chickens, skeletal muscle hyperplasia is thought to occur prior to hatch, and the increased size of the post-hatch muscle is largely due to hypertrophy [78]. Hypertrophy is driven by controlling the balance between protein synthesis and degradation [54,55]. The enrichment of 12 ribosomal subunit proteins, eukaryotic translation initiation factors EIF3I and EIF4A4, eukaryotic translation elongation factor 2EEF2, along with five tRNA ligases likely play an important role in the increased level of muscle hypertrophy seen in modern lines (Fig 2). Elevated protein synthesis could also cause the increase in stress response evidenced in the modern birds. The stress proteins function in a variety of processes including serving as chaperones or cochaperones, or in modulating genotoxic stress, redox regulation, and protein degradation (Fig 2).
The legacy line birds exhibited enrichment for HK1 and PFKM enzymes, which drive the first phase of glycolysis (Fig 3). PFKM catalyzes the rate limiting step of glycolysis, and its enrichment would direct Glucose 6-Phosphate down the glycolytic pathway. This is consistent with the typical fast-twitch breast muscle fibers that are seen in birds that exhibit brief episodes of flight. In contrast, modern birds are enriched in one protein (TKT) that supports the nonoxidative part of the pentose phosphate pathway (PPP). The non-oxidative portion of the PPP provides precursors for nucleotide synthesis and feeds glucose metabolites back into glycolysis. Modern Ross 708 birds also express higher levels of phosphoribosyl pyrophosphate synthetases, which direct Ribose-5-Phosphate to nucleotide production. The elevated level of LDHB in the Ross 708 birds is expected to drive the lactate:pyruvate equilibrium towards pyruvate. This would retain pyruvate for further metabolism by the TCA cycle and limit the release of lactate for energy production by other tissues [79]. Furthermore, glucose consumption supporting the modern broiler's breast muscle has ramifications for brain development. As glucose is the main energy source for the brain, diversion of this nutrient to the breast muscle may cause the reduced brain growth seen in modern broilers compared with legacy chicken lines or the red junglefowl [13].
The elevated levels of multiple enzymes of the TCA cycle in the Ross 708 breast muscle allows this pathway to meet demands of breast muscle hypertrophy (Fig 3). This is supported by studies in other species implicating elevated TCA cycle activity in hypertrophy. For example, enrichment of TCA metabolites was noted in a Klf10 mouse knockout model of soleus muscle hypertrophy [80] with similar results seen in aerobic exercise induced hypertrophy in humans [81,82]. Furthermore, resistance training induced hypertrophy in humans increases the activity of citrate synthase, the gateway enzyme to the TCA cycle [83].
Several enzymes involved in lipid beta-oxidation are also elevated in birds from the modern line (Fig 4]. Elevated expression of ACAD9 is particularly informative as this is the rate-limiting enzyme controlling lipid oxidation and elevated HADH activity also plays a role in skeletal muscle hypertrophy [83]. Increased ACAD9, HADH, and other enzymes involved in lipid oxidation indicates that modern Ross 708 birds are using lipid metabolism, in addition to the pentose phosphate and TCA pathways, to provide resources supporting expansion of the breast muscle.
Elevated lipid beta-oxidation in the breast muscle may have ramifications for morphometric changes seen in the growth of modern broilers. For example, the reduced normalized heart mass in the modern Ross 708 line compared with legacy UIUC birds could cause the cardiomyopathy seen in modern broilers. Also, if normalized spleen mass is viewed as a proxy for immune functions, the morphometric data indicates that immune function is significantly lower in modern lines compared with birds from the legacy line. Metabolically this may arise from the elevated lipid use in modern Ross 708 skeletal muscle. Cardiac muscle and the immune system use lipids as a major source of energy. Consequently, increased competition with skeletal muscle for lipids might inhibit heart growth and immune function seen in modern broilers. In addition to glucose the brain also readily uses ketone bodies, such as acetoacetate, to function. The elevated levels of ACAT2 in modern breast muscle may reduce the availability of ketone bodies for use by the brain.
These data support prior studies comparing the transcriptomes of 6-day old modern Ross 708 and legacy UIUC birds which concluded that the legacy line breast muscle was enriched for transcripts associated with glycolysis, while the transcriptome of the modern birds favored beta-oxidation [16]. Additionally, TKT transcripts were also elevated in a study of birds with high feed efficiency compared with low feed efficiency chickens [26]. Elevated levels of this enzyme seen in this proteome analysis provide further support for increased expression of these proteins improving feed efficiency.
Although these samples were obtained at day 6 post-hatch, there are already differences between the modern and legacy lines that have implications for the development of Wooden Breast Disease [28]. The elevated levels of stress proteins seen in breast muscle from the modern line provide compelling evidence that this tissue is undergoing a variety of stresses, likely due to its rapid growth. Oxidative stress is thought to be a major contributor to the development of this disease [28, 84,85] and five of the stress responsive proteins detected in this study, PIT54 protein (PIT54) [86,87], peroxiredoxins (PRDX1, PRDX4, PRDX6) and thioredoxin (TXN) play important roles in regulating oxidative stress. Gene Ontology analysis of our data also detected proteins associated with neutrophils in the mammalian immune system. Neutrophils are one of the earliest responders to inflammation and in birds the role of these phagocytic cells is filled by heterophils [77]. Histological examination of birds prone to Wooden Breast Myopathy revealed heterophilic infiltration in the pectoral muscle of D14 chickens and this is thought to be an early sign of disease development [25]. Furthermore, lipid metabolism has been shown to be altered in rapidly growing broilers that develop Wooden Breast Disease [88]. Taken together, the proteomic data suggests that Wooden Breast Disease starts to develop well before the disease is visually or palpably evident.

Conclusions
In support of our earlier transcriptome studies [16], these data indicate specific responses to artificial selection for improved feed efficiency and breast muscle yield in the modern broiler. The breast muscles of Ross 708 birds exhibit increased ribosomal protein synthesis, along with changes in the glycolytic, pentose phosphate, TCA, and beta-oxidation pathways when compared with the UIUC chickens. A reasonable hypothesis is that these adaptations are the direct result of artificial selection as they would support the improved feed efficiency and rapid breast muscle growth seen in the Ross 708 chickens. The selection possibly drove some of the negative outcomes seen in modern broilers. The diversion of nutrients to support the excessive breast muscle growth could cause the reduced growth of the brain and heart along with reduced immune function. In addition, the increased levels of proteins controlling oxidative damage and heterophil function suggest that the Wooden Breast myopathy seen at market age may begin very early post-hatch. Finally, the identity of these differentially regulated proteins generates testable hypotheses regarding the regulatory mechanisms that orchestrate the changes that have been introduced to the broiler chicken by human selection.
Supporting information S1 Table. Differentially expressed proteins from the legacy (UIUC) and modern (Ross 708) broiler lines. Included are proteins and gene accessions from NCBI and expression data supported by transcriptome expression [16] is indicated (RNASeq Concordance).   (1)